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Abstract 

Aims. We investigate the surface nitrogen content of the six magnetic O stars known to date as well as of the early B-type star r 
Sco. We compare these abundances to predictions of evolutionary models to isolate the effects of magnetic field on the transport of 
elements in stellar interiors. 

Methods. We conduct a quantitative spectroscopic analysis of the sample stars with state-of-the-art atmosphere models. We rely on 
high signal-to-noise ratio, high resolution optical spectra obtained with ESPADONS at CFHT and NARVAL at TBL. Atmosphere 
models and synthetic spectra are computed with the code CMFGEN. Values of N/H together with their uncertainties are determined 
and compared to predictions of evolutionary models. 

Results. We find that the magnetic stars can be divided into two groups: one with stars displaying no N enrichment (one object); and 
one with stars most likely showing extra N enrichment (5 objects). For one star (0 1 Ori C) no robust conclusion can be drawn due 
to its young age. The star with no N enrichment is the one with the weakest magnetic field, possibly of dynamo origin. It might be a 
star having experienced strong magnetic braking under the condition of solid body rotation, but its rotational velocity is still relatively 
large. The five stars with high N content were probably slow rotators on the zero age main sequence, but they have surface N/H 
typical of normal O stars, indicating that the presence of a (probably fossil) magnetic field leads to extra enrichment. These stars may 
have a strong differential rotation inducing shear mixing. Our results should be viewed as a basis on which new theoretical simulations 
can rely to better understand the effect of magnetism on the evolution of massive stars. 
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1. Introduction 

The evolution of massive stars is governed by two main physi- 
cal processes: their mass loss (Chiosi & Maeder 1986) and their 
rotation (Maeder & Meynet 2000). All additional processes di- 
rectly affecting these properties are thus major contributors to 
the fate of massive stars. This is the case of magnetism. Since the 
pioneering work of Babel & Montmerle ( 1 997), it is known that 
the presence of a large scale surface magnetic field will modify 
massive stars winds. By deflecting material along the field lines, 
magnetism affects the wind geometry and structure. Simulations 
by ud-Doula & Owocki (2002); Ud-Doula et al. (2008, 2009) 
have shown that equatorial overdensities can be created. Material 
can either fall back onto the stellar surface or be ejected in the 
magnetic equatorial plane depending on the interplay with ro- 
tation. In the case of very strong magnetic fields, observations 
and simulations show that a rigid magnetosphere strongly affects 
the wind structure, creating caps of stellar material in the min- 
ima of the magneto-rotational potential (Townsend et al. 2007). 
Since the wind flow is deflected by the presence of a magnetic 
field, the way the star loses its angular momentum is affected. In 
addition, the internal structure of the star is also affected, forc- 


* Based on observations collected at the CFHT and the Telescope 
Bernard Lyot. 


ing the star to follow new evolutionary paths. Given the rele- 
vance of massive star evolution to many fields of astrophysics 
(galactic chemical evolution, interstellar dynamics, supemovae 
and gamma-ray bursts, nucleosynthesis...), it is crucial to under- 
stand the magnetic properties of these objects in order to con- 
strain their influence on the evolution and fate of OB stars. 

The recent development of powerful spectropolarimeters 
has revolutionized our view of massive stars magnetism. The 
first detection of a magnetic field on the Trapezium O-fype 
star 0 1 Ori C (Donati et al. 2002) dates back to less than 
a decade. This spectroscopically varying and strong X-ray 
emitting star has revealed a surface kilogauss field organized 
mostly as a dipole. Subsequent detections have been made on 
HD 191612 (Donati etal. 2006a), HD 148937 (Hubrigetal. 
2008; Wade et al. 2011b), HD 57682 (Grunhut et al. 2009) and 
HD 108 (Martins et al. 2010). Field strengths are of several hun- 
dred Gauss to a few kilogauss and the geometry, when con- 
strained, is consistent with a dipole. The only exception is the 
O supergiant £ Ori A (Bouret et al. 2008) for which a more com- 
plex field topology was tentatively derived, together with a re- 
markably weak field of 50-100 G. This raises the question of 
the origin of magnetism in massive stars. The debate is ongoing. 
A fossil origin, stable over Myrs and with a rather simple ge- 
ometry, is supported by some simulations (Braith waite & Spruit 
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20G4; Duez et al. 2010). In that case, the star retains the orig- 
inal magnetic field of its parental molecular cloud or dynamo 
field of its convective progenitor. Alternatively, dynamo pro- 
cesses may be at work. In analogy with low mass stars where 
a ccnvective envelope (or even a fully convective structure) ex- 
ists, magnetic field might be generated in the convective core 
of high mass stars and subsequently transported to the surface 
(Charbonneau & MacGregor 2001; Bnin et al. 2005). A dynamo 
may also operate in the radiative zones of massive stars (Spruit 
2002; MacDonald & Mullan 2004). The field would be produced 
by the shear in this region. This idea has been tested in numer- 
ical simulations, with contradictory results (Braithwaite 2006; 
Zahn et al. 2007). 

However, assuming that this dynamo process was able to 
produce and maintain a magnetic field, evolutionary calcula- 
tions have been conducted to test the effects of this type of 
magnetism on the evolution of massive stars (Maeder & Meynet 
2004). They find that the internal structure is strongly affected, 
with a completely different rotation profile in the star compared 
to non-magnetic calculations including rotation. A consequence 
of this re-oiganization of the stellar interior is a modification of 
the transport of angular momentum and of chemical elements. 
In particular, Meynet & Maeder (2005) studied the effects on 
the surface chemical appearance of the star. They show that the 
presence of a dynamo generated magnetic field significantly in- 
creases the efficiency of mixing, leading to strong nitrogen and 
even helium surface enhancement. Further analysis of the feed- 
back effects of magnetic braking balanced these results though, 
highlighting the importance of the type of rotation profile (dif- 
ferential versus solid-body) in the final chemical appearance of 
the star (Meynet et al. 201 1). 

Although still limited, the sample of magnetic O stars is now 
populous enough to test these predictions. An important side ef- 
fect of the spectropolarimetric observations is the production of 
very high S/N ratio optical spectra at high resolution. These are 
perfectly suited for abundance studies. In this paper, we thus an- 
alyze the surface nitrogen content of the magnetic O stars known 
to date as well as of the B0.2V star r Sco. Our aim is to see if 
peculiar chemical patterns are observed and if so, to what extent 
they can be explained by the current generation of evolutionary 
models. The goal is to identify specific patterns that can only 
be attributed to the presence of a magnetic field. This is par- 
ticularly important in the context of explaining the population of 
N-rich slowly rotating stars in the Magellanic Cloud reported by 
Hun:er et al. (2008, 2009). Such stars are not explained by rotat- 
ing evolutionary tracks and have been suggested to be magnetic. 
The observations are presented in Sect. 2. We describe the analy- 
sis with atmosphere models in Sect. 3. The results are discussed 
in Sect. 4. We show how the chemical properties of the magnetic 
O stars can be interpreted in terms of evolutionary models with 
rotation and with or without magnetic field. We finally summa- 
rize our conclusions in Sect. 5. 


2. Stellar sample and observational data 

2.1. Observations 

The observations were collected with echelle spectropolarime- 
ters ESPaDOnS at CFHT (Canada-France Hawaii Telescope) 
and NARVAL at TBL (Telescope Bernard Lyot), partly from 
the MiMeS project (Wade et al. 201 la). The spectra cover the 
wavelength range 3700-10500 A with a resolution of 68000 
(ESPaDOnS) and 65000 (NARVAL). We gathered multiple se- 
quences for each of the sample stars, each sequence consisting 


of four individual subexposures taken in different polarimeter 
configurations. From these four individual subexposures, mean 
Stokes parameters were derived. All frames were processed us- 
ing Libre-ESpRIT (Donati et al. 1997), a fully automatic reduc- 
tion package installed at CFHT and TBL for optimal extrac- 
tion of ESPaDOnS and NARVAL spectra. The peak signal-to- 
noise ratios (SNRs) per 2.6 km s _1 velocity bin depends on the 
instrument, on the exposure time and on weather conditions. 
Complete references to the spectropolarimetric data are gathered 
in Table A. 1 . A summary of the data available is given in Tab. 1 . 
We used the average of the entire set of spectra for each star, en- 
suring that no variability in the key diagnostic lines was present 
among the data. 

2.2. Sample stars 

Our sample is composed of the six magnetic O stars known to 
date. Here we give a few details about each of these objects. The 
journal of observations is shown in Table A.l, 

HD 191612 is a member of the Of?p class, first intro- 
duced by Walbom (1972). Such stars exhibit, among other pe- 
culiarities, variability in their Balmer, He i, C m and Si m lines. 
Donati et al. (2006a) reported the detection of a kilogauss mag- 
netic field on its surface. It is presently well-established that 
HD 191612 displays two recurrent states, with a period of ap- 
proximately 538 days (Walbom et al. 2003; Walbom & Howarth 
2007; Howarth et al. 2007). These two states are referenced here 
as ’’quiet” (i.e. 08 spectral type, Ha in absorption) and ’’active” 
(i.e. 06 spectral type, strong wind signatures). We decided to 
focus on the quiet phase since the photospheric lines, used to 
derive the stellar parameters and surface abundances, are less 
contaminated by wind emission than in the active phase. We re- 
lied on data taken on August 23 rd 2009, corresponding to a phase 
of 0.393 (see Wade et al. 201 lc). The journal of observations in 
Table 1 presents additional information. 

HD 108 is the second Of?p star of our sample. It is very 
similar to HD 191612, with a long-term variability both photo- 
metrically and spectroscopically. However the timescale of this 
variability is much longer and is estimated to be of about 55 
years. The He i and Balmer lines are the most affected, leading 
to a spectral type changing from 04 to 08. Martins et al. (2010) 
reported the detection of a (most likely) kilogauss magnetic field 
on its surface. We used the average of the July- August 2009 sub- 
set of the Martins et al. observations in the present study since, 
as for HD 191612, this corresponds to the minimum of extra 
emission in the Balmer and He i lines. 

HD 148937 is the last member of the Of?p class in our sam- 
ple. Its variability is much more limited (lower amplitude) than 
that observed in HD 191612 and HD 108 spectra (Naz6 etal. 
2008b). The periodicity of the variations of about 7 days, is ex- 
tremely short compared to the long periods observed in the two 
other Of?p stars. The HD 148937 spectra were collected during 
six nights in June 2010. 

0* Ori C (HD 37022) is a mid-late type O dwarf. It is the 
first O star for which a magnetic field was detected (Donati et al. 
2002). It is known for its variability, with a period of approxi- 
mately 15.4 days (Stahl et al. 1996; Wade et al. 2006), the most 
striking variable spectral features being an inverted P-Cygni pro- 
file in the Hen/14686 line and strong emission features in the 
Balmer lines (especially Har and H0). The equivalent width of 
some He and metal lines also vary as a function of the phase 
(Simon-Diaz et al. 2006, and references therein). 

Spectropolarimetric observations of £ Ori A were collected 
with NARVAL in 2007 October. Altogether, 292 circular polar- 
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Table 1. Summary of observations. 


Star 

Instrument 

Total Exp. time 

Number of spectra 

typical S/N per spectrum 

HD 108 

NARVAL 

34 h 20 min 

37 

150-500 

HD 148937 

ESPADONS 

4 h 00 min 

6 

200-700 

HD 191612 

ESPADONS 

2 h 40 min 

8 

100-350 

O' OriC 

NARVAL 

53 min 

2 

500-1300 

£ Ori A 

NARVAL 

13 h 52 min 

336 

1000-1400 

HD 57682 

ESPADONS 

3 h 4 min 

5 

400-800 

r Sco 

ESPADONS 

13 min 

5 

600-1200 


Notes. The complete journal of observations is given in Appendix A. 1 


ization sequences were taken. The peak SNR per 2.6 km s _1 ve- 
locity bin varies between 800 and 1500. £ Ori A probably has a 
short rotation period of approximately 7 days and a weak mag- 
netic field (Bouret et al. 2008). Although some spectral lines ap- 
pear variable, the main diagnostics used to derive photospheric 
parameters remain stable during the rotation cycle. 

HD 57682 is an 09 sub-giant. Its magnetic field of strength 
1.7 kG was reported by Grunhutetal. (2009). HD 57682 dis- 
plays spectral variability but on a much lower level compared to 
HD 191612 and HD 108. Its rotational period is estimated to be 
~64 days (Grunhut et al. 201 1). In the present work, we used a 
set of exposures taken between May 4 th and May 9 th 2009. 

Finally, we have included r Sco (HD 149438) in our sample. 
Although not an O-type star per se, r Sco is an early B dwarf 
(spectral type B0.2V) just at die limit of the O star mass range. 
Since a magnetic field of the order 300 G has been clearly de- 
tected on its surface (Donati et al. 2006b), it was justified to add 
it tc our otherwise relatively limited sample. The geometry of 
the magnetic field is quite complex, with energy in multipoles 
up to order five. The origin of the field is debated, although a 
fossil origin is preferred. We have used a selection of spectra 
taken between 2005 and 2008 to conduct our analysis. 

To maximize the signal-to-noise ratios, we used the aver- 
age of the observations presented in Table A. 1 for all stars. This 
resulted in spectra with SNR above 1000. 

3. Spectral analysis 

3.1. Model Atmospheres 

The spectroscopic analysis was performed with the atmosphere 
code CMFGEN (a complete description of the code is provided 
by Hillier & Miller 1 998). The models include stellar winds and 
are computed in spherical geometry. The code solves simulta- 
neously the radiative transfer equation (in the co-moving frame) 
and the rate equations, following an iterative scheme. The com- 
putations provide non-LTE occupation numbers. The tempera- 
ture structure is given by the radiative equilibrium equation and 
the density structure by the mass conservation equation. 

Line-blanketing is included in CMFGEN through a super- 
level formalism in order to reduce the computational cost. Our 
models include H, He, C, N, O, Si, Ne, Mg, Si, S, Fe andNi. For 
each species, the solar abundances from Grevesse et al. (2007) 
were taken as reference. 

The wind structure is described by a mass loss rate, a 
terminal velocity and a /^-velocity law (Pauldrach et al. 1986; 
Kudritzki et al. 1989). The subsonic part of the atmosphere is 
described by a pseudo-hydrostatic structure, computed with the 
TLUSTY code (Lanz & Hubeny 2003) with a micro-turbulent 
velocity = 10 km s -1 . The velocity structure of the at- 


mosphere is then fully defined by connecting the subsonic 
and the supersonic components. We have used clumped mod- 
els, with a volume filling factor of 0.1 (0.01 for ( Ori - 
see Bouret et al. (2008) and HD108 - see Naze et al. (2008a)). 
The other wind parameters have been adopted from the litera- 
ture (Donati et al. 2006b; Bouret et al. 2008; Naz6 et al. 2008a; 
Grunhut et al. 2009; Wade et al. 201 lc,b). 

Once a model atmosphere is converged, the final synthetic 
spectrum is obtained from a formal solution of the radiative 
transfer equation. A depth-dependent micro-turbulent velocity 
£(r), equal to 5-10 km s" 1 at the base of the wind (f 0 ) and to 
about 10 % of the terminal velocity far from the photosphere 
(£»), is allowed. The value used for the models has very lit- 
tle effect on the atmosphere structure, and consequently on the 
emergent spectrum. On the other hand, fo may modify the shape 
and strength of the photospheric lines. We found that a value of 
1 0 km s _1 gave satisfactory fits for most stars, except HD 57682 
(the most sharp-lined star of the sample) for which a value of 
5 km s -1 was better. 

3.2. Stellar parameters 

The stellar parameters of our best models are gathered in Table 2. 
They correspond to the set of parameters that best fit the ob- 
served spectrum of the respective sample star. The uncertainties 
are indicated for each parameter. We will detail below how we 
determined the uncertainties on the nitrogen abundance. The er- 
rors on effective temperature and surface gravity depend on the 
quality of the fits we obtain, and are used to compute the errors 
on luminosity, radius and mass. The main source of uncertainty 
on the luminosity is the distance. 

We derived the effective temperature by the ionization bal- 
ance method. In practice, the relative strength of Hei and 
Hen lines was used. We relied preferentially on Hei ,14471 
and Hen ,14542. Additional diagnostics were: Hei /14026, 
Her ,14388, Hei 44712, Hei ,14920, Hei ,15876, Hen ,14200, 
He u ,15412. In the Of?p stars, the Hei lines are partly contami- 
nated by some additional emission the origin of which is unclear. 
The contamination is the largest in the “active” states, i.e. when 
the C III A 4650 lines have their maximum emission strength. 
Consequently, the temperature determination is difficult in those 
states and this explains why we rely on the quiet states whenever 
possible. 

The Balmer line wings are the best surface gravity diagnos- 
tics available in optical spectra. We relied mainly on H/J, Hy and 
H 8 to constrain logg. In the Of?p stars, the Balmer lines suf- 
fer from the same contamination as the He i lines in the active 
state. Fortunately, the wings, especially the blue one, are rela- 
tively unaffected by this contamination, allowing a safe gravity 
determination. Note that H 8 is blended with NIII ,14097, render- 
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Figure 1. Example of fit of the nitrogen lines in HD 57682. The observed spectrum is the black solid line. The smooth colored lines 
are models with the best fit parameters and different values of N/H. The^ 2 function is evaluated using these models to constrain the 
value of surface nitrogen content. In the present case, N/H = 1.3 x 10 -4 is the preferred value (see Fig. 2, red curve). More models 
with additional N/H values have been calculated but only four are shown for clarity. 


ing the use of this line less suitable. It is used as a secondary 
indicator. 

The luminosities are computed in two ways. If calibrated 
UV flux and optical- infrared photometry is available, the spec- 
tral energy distribution is fitted. The Galactic extinction laws of 
Seaton (1979) and Howarth (1983) are used to derive the extinc- 
tion The luminosity is adjusted to match the SED, for a given 
distance. If only optical photometry is available, we derived the 
luminosity from the V band magnitude, the effective tempera- 
ture. the bolometric correction (from Martins & Plez 2006), the 


extinction (calculated from B-V) and the distance. Once known, 
the luminosities directly provide the radii. 

The distance of the Orion stars is set to d = 414 ± 50 
pc (Mentenetal. 2007). The distance to the Of?p stars is de- 
termined by their supposed membership to OB associations, 
and was taken from Humphreys (1978). HD 108 belongs to 
the Cas OB5 association. We kept the luminosity derived by 
Martins et al. (2010) assuming a distance of 2.5 1±0. 10 kpc and 
fitting the spectral energy distribution. HD 148937 is a mem- 
ber the Ara OBI association, at a distance of 1.38±0.12 kpc. 
HD 191612 belongs to the Cyg OB3 association, at a distance of 
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2.29±0.12 kpc. r Sco has a reliable Hippacros parallax leading 
to a distance of 145± 1 1 pc. 

Grunhutet al. (2009) assumed the luminosity of HD 57682 
to be typical of a late O-type star and thus adopted the calibrated 
value of Martins et al. (2005). From this luminosity and the fit to 
the IUE spectrum, they derived a distance of 1 .3 kpc. We adopted 
their luminosity for our own models 

The projected rotational velocities and macroturbulent ve- 
locities have been determined from the fit of He i 4713. A first 
guess of V sin / was obtained from the line Fourier transform. It 
usually led to an upper limit, no zero being observed above the 
noise level. We then used various combinations ( V sinz / v mac ) 
to convolve synthetic spectra and selected the values best repre- 
senting the observed line profile (see also Howarth et al. 2007). 
We used a Gaussian profile to take macroturbulent velocities into 
account. 


3.3. Surface nitrogen abundance 

The main purpose of the present work is to derive surface nitro- 
gen abundances as accurately as possible, in order to compare 
them to the predictions from stellar evolution models. 

This task is challenging considering that the optical spec- 
tra of O stars exhibit a small number of nitrogen lines, and that 
most of these lines are relatively weak. In practice, we have used 
all nitrogen lines that our models were able to reproduce. Since 
we have a relatively heteregeneous sample of stars in terms of 
effective temperature, the number of lines used is different for 
each star. The only common lines are the Nni 4510 - 4525 fea- 
tures. These lines are relatively free of contamination and do not 
suffer from severe blending. Their formation do not seem to be 
affected by line-blanketing effects. Their lower levels are con- 
nected to the N in ground level by transitions at 434A, a region 
relatively free of metallic lines. Hence, subtle radiative transfer 
effects such as those reported by Najarro et al. (2006) are not ex- 
pected. For HD 57682 and r Sco we took advantage of the rela- 
tively large number of N n lines present in the observed spectra. 
For die Of?p stars, we only relied on N m features since the stars 
are too hot to display N n lines and N iv 4058 was clearly not 
reproduced. This inconsistency is probably not related to an ion- 
ization issue since for stars with lower T & the N m/Nn balance 
matches the Heu/Hei balances. Rivero Gonzalez et al. (2011b) 
reported that in some of their models, they could not correctly re- 
produce this feature (as well as other N rv lines). Part of the prob- 
lem might be due to the sensitivity of these lines to wind prop- 
erties (Rivero Gonzalez et al. 201 lb). We did not perform an in- 
vestigation of the wind properties of these objects, and thus we 
refrained from using this line when clear inconsistencies were 
observed. In addition, since magnetic field severely complicates 
the determination of wind properties (due to geometry, clump- 
ing), it is very likely that N rv 4058 is not perfectly predicted by 
our spherical models. Similar considerations can be applied to 
the strong N in 4634 - 4643 lines. Rivero Gonzalez et al. (2011a) 
showed that the onset of the wind was critical in the formation 
of this triplet. We note in addition that blends with O ii lines ex- 
ist a: 4643 A. We thus used these lines only when pure absorp- 
tion was observed. Besides, the Of?p stars are known to display 
narrow N m 4634 - 4643 emission lines of unknown origin (e.g. 
Naze et al. 2008a). In the case of 0 l Ori C, we could reliably use 
only the N m 45 1 0 - 4525 and N hi 4095 features (the latter being 
located in the wing of Hd and thus subject to more uncertainties). 
All the other N m features were too weak. 



N/H [x 10 4 ] 


Figure 2. Chi squares computed for the quiet state of HD 57682 
as a function of the nitrogen abundance. Each color represents 
a parameter set: red for the best fit model (parameters given in 
Tab. 2), green for variations of the surface gravity (± 0.1 dex 
with respect to the best fit model) and finally, blue for variations 
of the effective temperature (± 1 000 K with respect to the best fit 
model). The symbols correspond to the models computed. The 
different curves are polynomial fits to the respective set of sym- 
bols (red triangles, blue squares,. . .). The black dot with the error 
bars shows the position of the final N/H value and its uncer- 
tainty. 


Special care was taken to estimate the error bars on the de- 
rived surface nitrogen abundances. To do so, we proceeded as 
follows: for the best fit model (best values of and logg) we 
ran models at different nitrogen abundances. We then computed 
two new series of models with 7^ increased (respectively de- 
creased) by the uncertainty on 7^-. We proceeded similarly for 
two series of models with the best estimate of 7df, but logg in- 
creased (decreased) by the uncertainty on logg. We ended up 
with five series of models with 5-8 values of N/H for each se- 
ries. Fig. 1 shows an illustration of the model fits to the observa- 
tions in the case of HD 57682. 

We then estimated the preferred N/H content for each serie. 
For that, we calculated for the entire set of lines 


1 y M-Si? 

N ~ x it 


( 1 ) 


where M t is the model spectrum, S t the observed spectrum, cn 
the uncertainty on the observed spectrum and N the number of 
wavelength points over which the function is estimated. We de- 
termined the value of N/H that minimizes the chi-square for 
each set of parameters (T e n, log g). An illustration of x 2 versus 
N/H for the various models computed for HD 57682 is given in 
Fig. 2. 

To obtain the final value of N/H, we determined the best fit 
value for each combination (7^, logg) and computed the aver- 
age and one <j dispersion of this set of values. The final values 
are given in Table 2. 
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Table 2. Stellar parameters derived from quantitative analysis with atmosphere models. 


Star 

HD191612 

0 1 OriC 

; Ori A 

HD57682 

HD 108 

HD 148937 

T SCO 

Spectral type 

06-8f?p 

07V 

09.71b 

091V 

04— 8f?p 

05.5-6f?p 

B0.2V 

Ta (K) 

36000±2000 

38000±2000 

29500±1000 

34500±1000 

35000±2000 

40000±2000 

31000±1000 

logg 

3.75±0.10 

4.20±0.20 

3.25±0.10 

4.00±0.20 

3.50±0.10 

4.00±0.10 

4.10±0.10 

log(i/Lo) 

5.45±0.15 

4.80±0.13 

5.64±0.15 

4.79±0.25 

5.70±0. 10 

5.80±0.15 

4.37±0.08 

R (»©) 

12.3±1.9 

5.8±1.2 

25.4±1.7 

7.0±0.4 

19.4±2.2 

16.6±1.7 

5.3±0.6 

M»(M©> 

37.7£ 7 

23.8 + 2 9 

42.8J* 7 

22.0+ 3 / 

4S.81 4 

57 9* 7 0 

15.9!™ 

Mpec (Mq) 

38.7±20.7 

19.7±14.0 

42.0±21.2 

17.8±14.02 

43.3±18.7 

101.0±53.1 

13.1 ±4.5 

V sin; (km s -; ) 

< 1 

24 

100 

15 

< 1 

<50 

6 

v,,*,. (km s‘ l ) 

40 

30 

~0 

16 

45 

60 

3 

N/H (number) 

(2.7±L5) 10" 4 

(5.1 ±4.5) 10“ 5 

(3.3±1.0) 10" 5 

(1.3±0.3) 10" 4 

(5.9±4.2) 10" 4 

(3.0±1.2) 10” 4 

(l.4±0.2) 10“ 4 

12+log (N/H) 

8.43±0.24 

7.82±0.29 

7.52±0.14 

8. 1 1 ±0. 10 

8.77±0.31 

8.48 ±0.17 

8.15±0.06 


Notes. The evolutionary masses are determined using the tracks of Brott et al. (2011) with an initial rotational velocity of 300 km s l . 


To further check the robustness of our determinations, we 
ran a few additional tests. In particular, wc tested the effect of a 
change of the microturbulent velocity in the atmosphere model 
computation (not just the output synthetic spectrum). Increasing 
this velocity from 5 to 10 km s _l translates to a decrease of the 
derived abundance by 15-20%. Although not negligible, this is 
within the uncertainties of our determinations. 

Abundance determination have been performed by Morel 
(2011) for two our our sample stars: HD 57682 and t Sco. 
For die latter, our results are in very good agreement 
(12+log(W///)=8.15±0.06 versus 8.15±0.20). Przybilla et al. 
(2010) also found a value of 8. 16±0. 12, in good agreement with 
our findings. For HD 57682, we derive a nitrogen content much 
larger than Morel (201 1) (8.1 1 ±0.10 versus 7.52±0.25). We at- 
tribute partially this discrepancy to the use of Kurucz models 
in the Morel et al. study. They did not self-consistently derive 
the stellar parameters and the abundances, but adopted and 
logg from the study of Grunhut et al. (2009). Another cause of 
the problem can be a different line list used to perform the anal- 
ysis. However, we do not have access to the lines used by Morel 
et al. and cannot test this possibility. 

4. Stellar magnetism and nitrogen enrichment 

4.1. Comparison to predictions of rotating models without 
magnetic field 

In Fig. 3 we show the Hertzsprung-Russel (HR) diagram with 
the positions of the magnetic OB stars analyzed in the present 
study. We have chosen the evolutionary tracks and isochrones 
of Brott et al. (201 1). Other evolutionary tracks exist in the lit- 
erature, but the Brott et al. tracks have been computed for a 
large set of initial rotational velocities, allowing an investiga- 
tion of the effects of rotation on the results. In addition, they 
have the same initial metal content as that used in our atmo- 
sphere models, rendering possible comparison between our re- 
sults and the predictions of evolutionary calculations. Compared 
to the recent Geneva tracks (Ekstrom et al. 2011) they predict 
a stronger nitrogen enrichment at a given age (for moderate and 
zero rotation) during the early evolution representing our sample 
stars. Consequently, our conclusions regarding a possible extra 
enrichment of magnetic stars would be rc-enforced if we were 
using the Geneva tracks (see below). In Fig. 3, we have selected 
the evolutionary tracks with an initial rotational velocity of 300 
kms" 1 . This value is typical since it is supposed to lead to an 
average rotation on the main sequence similar to that observed 
(Maeder & Meynet 2000). 



Figure 3. Hertzsprung-Russel diagram for the sample magnetic 
OB stars. The isochrones and the evolutionary tracks are from 
Brott et al. (201 1) and were computed with an initial rotational 
velocity of 300 km s" 1 . 


From the HR diagram, one can estimate the evolutionary 
mass of each individual star. This is done by simple interpolation 
between the available tracks. The results are given in Table 2. 
The uncertainties reflect the errors on the effective temperature 
and luminosity. The age of the stars can also be tentatively de- 
rived. The ages are just a way of evaluating the degree of evolu- 
tion of a given star. They are directly used to build Fig. 4 where 
we show the surface nitrogen content as a function of age. This 
plot, together with the HR diagram, can tell us if a star is evolv- 
ing according to the predictions of standard evolutionary models, 
or if significant deviations are detected. 

Let us take HD 148937 as an example. From the HR dia- 
gram, its evolutionary mass should be 57.9 ^jq 4 Mq. The star is 
2.8^ £ Myr old. Looking at Fig. 4, we see that for this range of 
ages, the surface nitrogen content of HD 148937 is consistent 
with stars of initial mass 60 Mq. This is thus consistent with the 
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Table 3. Lines used to perform Nitrogen abundance determination. 



HD 191612 

0‘ OriC 

£ Ori A 

HD57682 

HD 108 

HD148937 

r Sco 

Nil 3995 


- 

~^r~ 

~T~ 

- 


V 

Nra 4004 

- 

- 

~ 

~ 

— 

V 

V 

Nil 4035 

— 

* 


V 
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Nil 4041-4043 

— 

- 

V 

V 



V 

Niv 4058 

• 

— 


V 

- 

- 

V 

Nm 4095 

— 

V 

- 
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— 


N in 4196 

V 


V 

V 

V 


V 

N m 4216 

V 

— 


V 

V 

- 

V 

Nn 4447 

— 

- 

V 

V 



V 

N m 4510-4525 

V 

V 

V 

V 

V 

V 

V 

N hi 4537 

V 



- 

V 

V 

V 

Nil 4602 

- 

- 

— 


* 


V 

Nu 4607 

- 

— 

— 

V 

- 

— 

V 

Nn 4613 


• 

— 

V 

— 

— 

V 

Nil 4621 

— 

- 

- 

V 


•»- 

V 

Nn 4634-4643 


- 

— 

V 

— 

— 

V 

Nn 4788 

— 

• 

- 

V 

• 

— 

V 

Nil 4803 

— 

- 

— 

V 


— 

V 

N hi 4907 


• 

— 

V 

— 

— 


Nn 4995 

- 

- 

- 

V 

- 

_ 

V 

Nn 5000-5011 



V 

V 

— 

— 

V 

N n 5026 

— 

- 




— 

V 

N n 5045 

— 

- 

— 

V 

• 


V 

Nn 5676-5679 



V 

V 

- 

- 

V 



Figure 4. N/H as a function of age. Black solid lines are the pre- 
dictions of evolutionary models including rotational mixing with 
an initial rotation rate of 300km s' 1 . Models are from Brott et al. 
(2011). The symbols are the magnetic O stars studied in the 
present paper. 


estimated initial mass of HD 148937, from which one can con- 
clude that this star does seem to have normal surface abundance 
for its evolutionary status. A similar conclusion can be drawn for 
HD 191612 and HD 108, although for the latter a hint of extra 
enrichment is not excluded. 


r Sco has an initial mass of about 16 M©. Compared to the 
15 M© evolutionary tracks, its nitrogen content is rather high. 
Only for an effective temperature at the low limit of the derived 
range (and thus a relatively old age) can the N content be barely 
reconciled with the predictions. HD 57682 also has a nitrogen 
content 2 to 3 times the initial value. This value is marginally 
compatible with a 22 M© star provided its age is on the upper 
side of the range of possible values. Hence, r Sco and HD 57682 
quite possibly show some extra enrichment. 0 1 Ori C is very sim- 
ilar in mass but it is too young to have produced a significant 
amount of nitrogen and thus the low N/H we derive is not unex- 
pected. 

The last star of our sample, £ Ori A, is the oldest and most 
evolved of the magnetic O stars. At an age of 3.7±0.6 Myr and 
with an initial mass of 42.8^*][, its surface N/H should reach at 
least 2.0 x 10" 4 . We derive a value about 4 times lower (4.4 ± 
1.0 x 10“ 5 ). As already noted by Bouret et al. (2008), ( Ori A is 
thus particularly chemically unevolved for its age. 

So far, we have compared the properties of the magnetic O 
stars to evolutionary models with an initial rotational velocity 
of 300 km s -1 . This value leads to an average velocity on the 
main sequence consistent with observations of normal O stars. 
However, one might wonder if this value is appropriate for mag- 
netic stars. The presence of a magnetic field is known to reduce 
the rotation rate as illustrated by Townsend et al. (2010) in the 
case of cr Ori E. Donati et al. (2006a) and Martins et al. (2010) 
have also noted the extremely low rotation of HD 191612 and 
HD 108. The question is: were these stars slow rotators already 
on the zero age main sequence? If so, the comparison with low 
V sin / evolutionary tracks is necessary. Table 4 presents the spin 
down timescales due to magnetic braking calculated according 
to Ud-Doula et al. (2009). We have also listed the ages esti- 
mated from Figs. 3 and 5. As noted by Donati et al. (2006a) 
and Martins etal. (2010), HD 191612 and HD 108 have spin 
down timescales marginally compatible with their age, implying 
that they might have experienced a reduction of their rotation 
rate by a factor 2-3. However, given their very low present-day 
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Table 5. Evolutionary masses derived for different initial rota- 
tional velocities. 


Star 

M„ (300) 

(100) 

HD 191612 

37-7 

57.9*?.. 

37 5* 61 

HD 148937 

HD 108 

48.8 + j« 

47. 

©' OriC 

23.8-lf 

23.6! : 

( Ori A 

42 - 8 4 : 

22.0+f 
15 9 + «’ 6 

43.0^' 

21.5+f 6 

HD 57682 

rSeo 

15.7*7 


Notes. The models used are from Brott et al. (201 1). 


V sin/, those stars probably had initial projected rotational veloc- 
ities not higher than a few tens of km s _1 . All the other magnetic 
O stars have spin down timescales longer than their age. Given 
their rather low V sin / (with the exception of ( Ori A), these 
stars were also most likely slow rotators on the zero age main 
sequence. The youngest star of our sample (0'Ori C) confirms 
this trend: it is located close to the ZAMS and features a low 

V sin / (by O star standards). 

In view of these comments, we show in Fig. 5 the HR and 
N/H-agc diagrams with evolutionary tracks computed for an 
initial rotational velocity of 100 km s" 1 . Using these plots, we 
can see how the effects of rotation affect the previous results. 
The evolutionary masses derived from these new tracks are given 
in Table 5. The three Of?p stars (HD 108, HD 148937 and 
HD 191612) all appear to be much more N-rich than predicted 
by the 60 M© evolutionary track, and thus more enriched than 
expected for their mass, r Sco is now very clearly enriched com- 
pared to the 15 M© tracks which barely departs from the initial 
N content in the first 6 Myr. For HD 57682, an extra enrichment 
is also clearly present. For Q l Ori C, again the young age pre- 
vents any firm conclusion other than that no dramatic nitrogen 
enrichment is observed. For £ Ori A the surface chemical pat- 
tern is lower than predicted by the evolutionary tracks although 
marginal agreement exists if extreme stellar parameters are used. 

Given the above discussion, the following conclusions can 
be drawn: 

C Ori A: For moderate rotation rate, the low N content cannot 
be explained by standard mixing. Some extra process lead- 
ing to the quenching of rotational mixing is necessary to ex- 
plain the peculiar N/H ratio. The conclusion remains valid 
for low rotation rates, although the agreement with the ex- 
pected N/H ratio is marginal given the uncertainties on the 
age and surface nitrogen content. 

- HD 108, HD 148937, HD 191612, r Sco and HD 57682: 
under the reasonable assumption that these stars were slow 
rotators on the zero main sequence (i.e. initial rotation ve- 
locity of a few tens of km s"D, all stars show a higher de- 
gree of nitrogen enrichment than non-magnetic evolutionary 
models predict. In the less likely case where the initial rota- 
tional velocity was normal (i.e. around 300 km s~*), HD 1 08, 
HD 57682, r Sco and HD 191612 show marginal evidence 
for extra enrichment, while HD 148937 appears normal. 

— O^ri C: this star does not show peculiar nitrogen enrich- 
ment whatever the initial rotational speed. This is most likely 
because it is too young. 


4.2. Comparison to non-magnetic stars 

To further investigate the surface chemical properties of the six 
magnetic O stars, we show in Fig. 6 the same HR diagram as in 
Fig. 3 to which we have added comparison O stars. The latter 
objects have been taken from Bouret et al. (2005), Martins et al. 
(201 1), Bouret et al. (2012, A&A, in prep), Escolano et al. (2012, 
A&A, in prep), Przybilla et al. (2010) and Hunter et al. (2009). 
All O stars have been analyzed with CMFGEN in the same way 
as in the present study. In Fig. 6 we encode the value of N/H 
by the size of the symbol for each star (including the magnetic 
O stars). In the upper part of the HR diagram, HD 108 and 
HD 148937 fall in a densely populated group of O supergiants 
with N/H larger than 2.0 x 10 4 . They thus appear as similarly 
enriched compared to these objects, and present no peculiarity 
in their chemical pattern provided they were initially rotating at 
about 300 km s _1 (but see below). 

HD 191612 is located at a lower luminosity than the two 
other Of?p stars. In that region of the HR diagram, the surface 
nitrogen abundance of the few O stars present a large range of 
values, from less than 10 -4 to more than 4.0 x 10 -4 . In addition, 
the less evolved stars are more enriched than the more evolved 
objects. Given these puzzling results, the surface abundance of 
HD 191612 do not appear to deviate from the general trend, but 
the large range of N/H values observed in that part of the dia- 
gram prevents any firm conclusion regarding its ’normality’. 

The last star in the upper part of the HRD is £ Ori A. Its 
surface nitrogen content strongly deviates from the values ob- 
served in neighboring stars. The O supergiants analyzed so far 
all show significant N enrichment (N/H > 2,0 x 10~ 4 ) while £ 
Ori A is basically chemically unevolved. This confirms our pre- 
vious conclusions that this star clearly shows peculiar chemical 
properties. 

Going down in the HR diagram, we find 0‘Ori C and 
HD 57682 very close to each other. Only three O stars are lo- 
cated close enough to allow a quantitative comparison. These 
objects show a wide range of N/H values, from less than 1.0 
1CT 4 to more than 4.0 x 1(T 4 . Both O 1 Ori C and HD 57682 show 
relatively modest enrichment This is consistent with the two less 
enriched comparison stars. Hence, the two magnetic stars do not 
seem to stand out as peculiar objects. Finally, close to the loca- 
tion of r Sco we have found very few comparison stars and so 
we cannot draw any relevant conclusion. 

At first sight, this comparison only exacerbates the peculiar 
properties of £ Ori A. All the other stars appear similar to non- 
magnetic O stars. There is however an important property of the 
magnetic stars to consider: all (except £ Ori A) are slow rotators. 
As explained above, there are indications that they might have 
been slow rotators already on the ZAMS. As such, they are dif- 
ferent from the comparison stars which have diverse values of 
V sin / (32 stars with V sin / from roughly 0 to 330 km s - * 1 with 
an average of 116 km s _l ). Consequently, finding nitrogen en- 
richment in HD 108 and HD 148937 (the two objects for which 
there are sufficient comparison stars to draw firm conclusions) 
similar to that of comparison stars while they are rotating more 
slowly can be interpretated as the sign of an especially efficient 
mixing inside the star. Said differently, the (presumably) initially 
slowly rotating magnetic stars have the same surface N content 
as normal O stars. This is an indirect indication that the presence 
of a magnetic field in those objects tends to favor the transport 
of nitrogen from the core towards the surface. 

Clearly, a larger sample of magnetic and comparison stars 
would be needed to strengthen these conclusions. But the above 
comparisons all indicate that there seems to be two groups of 
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Table 4. Magnetic spin-down timescales and ages. 


Star 

Bp 

[kG] 

M 

M 0 

R 

R© 

Voo 

km s~ ] 

M 

[10“ 7 Mq yr -1 ] 

r’spin 

[Myr] 

age(300) 

[Myr] 

age(100) 

[Myr] 

HD 191612 

2.5 

37.7 

12.3 

2400 

0.3 

3.8 

2.8 + “ 

2.6+“ 

HD 148937 

1.0 

57.9 

16.6 

2600 

0.3 

11.3 

2.1+fl 

2.0*H 

HD 108 

0.5-2 

48.8 

19.4 

2000 

1.0 

2.0-7. 8 

2.9f 

2.0! : 

0 1 OriC 

1.1 

23.8 

6.5 

2500 

3.0 

3.24 

o.o+f! 

0.0+?! 

£ Ori A 

0.05-0.1 

42.8 

25.4 

2100 

17.0 

6.5-13.0 

3.6! : 

3.6! 

HD 57682 

1.7 

22.0 

7.0 

1200 

0.014 

18.8 

2.3+f 

2.8*|| 

r Sco 

0.5 

15.9 

5.3 

2000 

0.3 

17.0 

[ 3+33 
. 1:5 -2.3 

2 1 + ” 
Za -2.6_ 


Notes. The magnetic spin-down timescale is calculated according to Ud-Doula et al. (2009). Bp is the polar magnetic field. The stellar and wind 
parameters are used to calculate T sp i n . We have assumed a value of 0.1 for the moment of inertia constant k (see Eq. 25 of Ud-Doula et al. 2009). 
The entry ’age(300)’ is the age of the star estimated from the HR diagram built with evolutionary tracks computed for an initial rotational velocity 
of 300 km s -1 . Similarly, ’age(lOO)’ is the age for V sin* = 100 km s" 1 




(a) v=100 kins" 1 (b) v=100 km s" 1 

Figures. Same as Figs. 3 and 4 but for initial rotational velocities of 100 km s -1 . The initial masses derived from the evolutionary 
tracks are gathered in Table 5. 


magnetic stars: on one side £ Ori A which is clearly peculiar, 
with no nitrogen enrichment in spite of its evolved status and 
relatively normal rotational velocity; and the other OB stars, es- 
pecially the Of?p objects, possibly displaying larger enrichment 
thannon-magnetic stars provided they were indeed slow rotators 
ontheZAMS. 

4,3. On the nature of magnetic field in O stars 

The origin of magnetic field in massive stars is a matter of de- 
bate. The two main possibilities are either a fossil field or a 
field generated by dynamo processes in the radiative region or 
the convective core. In the former case, the field might come 
from the parental molecular cloud, the convective phases during 
pre-main sequence evolution, of from binarity. Relatively sim- 
ple geometries are expected in that case. For dynamo generated 
magnetic fields, several scenarios exists but the Tayler-Spruit 
dynamo has been the most tested (Spruit 1999, 2002). In this 
case, an initial poloidal field is transfered into a toroidal field 
by differential rotation. Later on, the Pitts-Tayler instability de- 


velops in the toroidal field and gives birth to a poloidal field. 
Consequently, the dynamo loop is closed. Controversy exists at 
present to know if this process is effectively at work. In a series 
of magneto-hydrodynamicai (MHD) simulations, Braithwaite 
(2006) claim to have observed the onset of the Pitts-Tayler in- 
stability and the dynamo loop. But Zahn et al. (2007) conclude 
differently in their simulations: although the Pitts-Tayler insta- 
bility seems present, it has no effect on the poloidal field which 
slowly decreases due to Ohmic decay. Hence the dynamo loop 
is not closed. A good summary of these processes is provided by 
Maeder (2009) and Walder et al. (2011) (see also Zahn 2011). 

Can we provide any constraint on the nature of O star mag- 
netism? Meynet & Maeder (2005) ran evolutionary calculations 
under the assumption that the Tayler-Spruit dynamo was work- 
ing. In that case, they could show that one of the main effects 
was an increase of the mixing and of the efficiency of chemi- 
cal element transport. The reason is the nearly solid body rota- 
tion achieved in the radiative zone in the presence of magnetic 
field. In non magnetic rotating models, chemical elements are 
mainly transported by shear turbulence caused by differential 
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Figure 6. Same as Fig. 3 with comparison stars added (grey 
squares). The comparison stars are from Bouret et al. (2005), 
Martins etal. (2011), Bouret et al. (2012, A&A, in prep), 
Escolano et al. (2012, A&A, in prep), Hunter et al. (2009) and 
Przybilla et al. (2010). The size of the symbols (including the 
six magnetic O stars) is now proportional to the N/H value. The 
evolutionary tracks are for v = 300 km s“ l . 

rotation. Meridional circulation - created by thermal imbalance 
due to gravity darkening - contributes little to the transport of 
elements, but is the main way of redistributing angular momen- 
tum. If the star rotates nearly solidly, shear turbulence is reduced, 
but meridional circulation is greatly enhanced so that the net ef- 
fect on the transport of chemical species is a much stronger mix- 
ing. Meynet & Maeder (2005) report an increase of the surface 
N/H of a factor 2 to 5 larger than in non magnetic rotating mod- 
els for a 15 Mo star of age 5 to 10 Myr. Maeder (2009) note 
that these results depend on the details of the dynamo process 
at w ork and on the strength of the magnetic field. However, the 
trend seems to be an significant over-enrichment compared to 
standard rotating models. In their more recent investigation of 
the role of magnetic braking, Meynet etal. (2011) note that if 
the star is rigidly rotating, braking is so strong that it suppresses 
all surface chemical enrichment. Thus the exact role of magnetic 
field on surface abundances is not clear. 

Our results are meant to provide observational constraints 
that could trigger further theoretical developments. At present, 
they can confront theoretical predictions as follows: 

- dynamo fields: the evolutionary calculations focus on the 
Tayler-Spruit dynamo as explained above. One of the ex- 
pectations of dynamo fields is that, similarly to the sun, they 
show complex morphologies and not simple dipolar struc- 
tures. The constraints on the field geometry of magnetic O 
stars are still scarce, but f Ori A stands out one more time 
(Bouret et al. 2008) since the first observations indicates the 
possible presence of several magnetic spots. If a dynamo 
process is producing the magnetic field of that star, then 
our results indicate that no surface enrichment is observed. 
According to Meynet etal. (2011) - (see also Zahn 2011) 


- this could be the sign of a nearly solid-body rotation and 
strong braking (see above). However, the projected rotational 
velocity of f Ori A is not especially low (V sin i = 1 10 km 
s" 1 ). This is not consistent with a strong braking unless the 
star was initially rotating very fast. But in that case, would 
the high level of rotation be compatible with the absence of 
surface chemical enrichment? Clearly, dedicated simulations 
are necessary to see if a star with a dynamo field can still be 
rotating relatively normally and not have produced any sur- 
face nitrogen excess after several Myrs. 

- fossil fields: it is usually thought that fossil fields have 
simple dipolar geometries (Donati & Landstreet 2009). 
Observations indicate that O^ri C, HD 57682, HD 191612 
and HD 148937 have dipolar fields. For HD 108, the phase 
coverage is too poor to provide any constraint, r Sco fea- 
tures a more complex magnetic geometry, but Donati et al. 
(2006b) attributes this to the youth of the star. If a dipole 
should remain after a few Myrs in fossil field magnetic 
stars, higher orders are expected to remain in the early 
evolution, after the formation from the parental molecu- 
lar cloud. Donati et al. (2006b) thus favoured a fossil ori- 
gin for the magnetic field of r Sco. According to Sect. 4, 
all those objects were initially slow rotators on the ZAMS. 
Consequently, the observed surface abundances arc higher 
than they would be if there was no magnetic field. In prac- 
tice, the degree of enrichment corresponds to that observed 
in stars with initial V sin i of 300 km s“\ while the initial 
rotational velocity of the magnetic stars was more likely of a 
few tens of km s" 1 . Unfortunately, there are no simulations 
of the effect of fossil fields on the chemical appearance of 
massive stars to which we can compare these results. We can 
thus simply state that dipolar, probably fossil fields tend to 
lead to extra mixing in the interior of O stars. According to 
the results of Meynet et al. (201 1), this requires the presence 
of strong shear created by differential rotation. Magnetic O 
stars with fossil fields might thus have a steep velocity struc- 
ture. 

Clearly, the effects of magnetic field on the properties of O 
stars are still poorly understood. A larger sample of magnetic 
stars would be needed to see if the apparent separation into two 
groups (slowly rotating enriched stars / normally rotating non- 
enriched stars) is confirmed. Analysis of more comparison stars 
is required too. In parallel, new MHD simulations and evolu- 
tionary models, especially including fossil fields, are necessary 
to better understand the physics at work in the interior and atmo- 
sphere of magnetic O stars. 

5. Summary and conclusions 

We have conducted a quantitative spectroscopic analysis of the 
six magnetic O stars known to date together with the early B 
star r Sco. Spectra collected with ESPADONS at CFHT and 
NARVAL at TBL have been used. Atmosphere models com- 
puted with the code CMFGEN have been computed. We have 
determined the stellar parameters and the surface nitrogen abun- 
dance. A careful investigation of the uncertainties associated 
with N/H has been pursued. 

Our results can be divided into two parts: 

• The surface N/H of £ Ori A - the star with the weakest 
magnetic field - is lower than the predictions of evolution- 
ary models. For that star the presence of the magnetic field 
probably quenches the transport of chemical elements. The 
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conditions under which such a process can happen are not 
clear. Evolutionary models accounting for magnetic braking 
predict no chemical enrichment in solid body rotating stars. 
But their rotational velocities are also rather low, whereas 
f Ori A appears to have a rather standard V sin i, 

• The other magnetic O stars, especially the three Of?p stars, 
display surface nitrogen content consistent with those of 
non-magnetic O stars with initial V sini of 300 km s -1 . 
However, these magnetic O stars were most likely slow rota- 
tor on the zero age main sequence (with values of V sini of a 
few tens of km s' 1 at most). Consequently, they feature larger 
surface N/H than they should if they did not host magnetic 
fields as also indicated by their nitrogen excess compared to 
rotating evolutionary models. 

The presence of a magnetic field in O stars thus appears 
to produce two types of chemical peculiarities: either an extra 
enrichment, or the absence of enrichment. Confrontations with 
evolutionary predictions in the case of dynamo fields indicate 
that solid body rotation might explain the properties of the non- 
enriched star, but its relatively normal rotational velocity is a 
puzzle in that context. For the N-rich magnetic O stars, a fos- 
sil field is usually assumed. Our results should be a basis for 
comparison with future evolutionary calculations including fos- 
sil fields since such evolutionary models do not exist at present. 
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Appendix A: Journal of observations 

Table A. 1 provides the complete journal of observations. 
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Table A.l. Journal of observations. 


Star 

Instrument 

date 

UT 

Exp. time (sec) 

S/N 

HD 108 

NARVAL 

Jul 21" 2009 

00:52:37 

4x750 

150-350 


NARVAL 

Jul 21 s1 2009 

02:18:22 

4x750 

200-500 


NARVAL 

Jul 2 1 st 2009 

03:11:49 

4x750 

200-500 


NARVAL 

Jul 25 th 2009 

00:41:56 

4x750 

200-500 


NARVAL 

Jul 25 th 2009 

01:35:23 

4x750 

200-500 


NARVAL 

Jul 25 lh 2009 

02:28:49 

4x750 

200-500 


NARVAL 

Jul 26 th 2009 

01:47:60 

4x700 

100-350 


NARVAL 

Jul 26 th 2009 

02:38:05 

4x700 

100-400 


NARVAL 

Jul 26 th 2009 

03:28:11 

4X700 

150-450 


NARVAL 

Jul 27 th 2009 

00:53:23 

4x806 

100-300 


NARVAL 

Jul 27 th 2009 

01:50:33 

4x806 

120-370 


NARVAL 

Jul 27 th 2009 

02:47:43 

4x806 

100-350 


NARVAL 

Jul 28 th 2009 

00:51:16 

4x900 

100-300 


NARVAL 

Jul 28 th 2009 

01:54:42 

4x900 

100-400 


NARVAL 

Jul 28 th 2009 

02:58:06 

4x900 

150-450 


NARVAL 

Jul 29 th 2009 

00:39:56 

4x900 

200-550 


NARVAL 

Jul 29 th 2009 

01:43:21 

4x900 

2504)00 


NARVAL 

Jul 29 th 2009 

02:46:46 

4x900 

200-500 


NARVAL 

Jul 30 th 2009 

00:56:23 

4x900 

250-630 


NARVAL 

Jul 30 th 2009 

01:59:50 

4X900 

250-600 


NARVAL 

Jul 30 th 2009 

03:03:16 

4x900 

250-600 


NARVAL 

Jul 31 s ' 2009 

00:30:28 

4x812 

200-500 


NARVAL 

Jul 31 st 2009 

01:28:02 

4x812 

200-500 


NARVAL 

Jul 31 st 2009 

02:25:37 

4X812 

200-550 


NARVAL 

Jul 31 st 2009 

03:23:11 

4x812 

200-500 


NARVAL 

Aug 1“ 2009 

00:55:57 

4x812 

200-550 


NARVAL 

Aug l a 2009 

01:57:06 

4x812 

250-600 


NARVAL 

Aug I s1 2009 

02:54:39 

4 X 812 

200-470 


NARVAL 

Aug 2 nd 2009 

00:50:14 

4x900 

200-470 


NARVAL 

Aug 2 nd 2009 

01:53:47 

4x900 

150-400 


NARVAL 

Aug 2 nd 2009 

02:57:24 

4x900 

150-400 


NARVAL 

Aug 4 th 2009 

01:24:42 

4x900 

200-470 


NARVAL 

Aug 4 th 2009 

02:28:07 

4x900 

250-400 


NARVAL 

Aug 4 th 2009 

03:31:34 

4x900 

200-570 


NARVAL 

Aug 5 th 2009 

00:22:43 

4x900 

200-500 


NARVAL 

Aug 5 111 2009 

01:26:09 

4x900 

200-600 


NARVAL 

Aug 5 th 2009 

02:29:35 

4x900 

200-600 

HD 148937 

ESPADONS 

Jun 20“ 2010 

09:10:53 

4x600 

300-850 


ESPADONS 

Jun 21 st 2010 

08:31:20 

4x600 

250-750 


ESPADONS 

Jun 22 nd 2010 

09: 9:60 

4x 600 

170-550 


ESPADONS 

Jun 23” 1 2010 

08:11:34 

4x600 

150-400 


ESPADONS 

Jun 24 th 2010 

08: 8:21 

4x600 

300-800 


ESPADONS 

Jun 25 th 2010 

08:16:39 

4x600 

3250-800 

HD 191612 

ESPADONS 

Aug 23“ 2008 

09:43:36 

4x300 

100-300 


ESPADONS 

Aug 23 rd 2008 

10:04:20 

4x300 

100-350 


ESPADONS 

Aug 23 ri 2008 

10:25:05 

4x300 

100-350 


ESPADONS 

Aug 23 ri 2008 

10:45:49 

4x300 

100-350 


ESPADONS 

Aug 23 ri 2008 

11:07:00 

4x300 

100-350 


ESPADONS 

Aug 23 rd 2008 

11:27:45 

4x300 

100-350 


ESPADONS 

Aug 23" 1 2008 

11:48:29 

4 X 300 

100-350 


ESPADONS 

Aug 23 rd 2008 

12:09:13 

4x300 

100-350 
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Table A.l. Continued. 


Star 

Instrument 

date 

UT 

Exp. time (sec) 

S/N 

0 1 OriC 

NARVAL 

Mar 10“ 2007 

08:28:19 

4x400 

500-1400 


NARVAL 

Mar 10* 2007 

08:59:01 

4x400 

450-1250 

£ Ori A 

NARVAL 

Oct 18'" 2007 

00:52:56 

48 x 4 x 20 

780-990 


NARVAL 

Oct 19 ,h 2007 

04:35:04 

8 x 4 x 40 

1010-1080 


NARVAL 

Oct 20 th 2007 

00:59:39 

44 x 4 40 

1220-1470 


NARVAL 

Oct 21" 2007 

23:06:15 

48 x 4 x 40 

810-1460 


NARVAL 

Oct 22 nd 2007 

23:45:52 

48 x 4 x 40 

1090-1480 


NARVAL 

Oct 23" 1 2007 

23:13:45 

48 X 4 x 40 

1030-1480 


NARVAL 

Oct 24 lh 2007 

23:58:34 

48 x 4 x 40 

1200 - 1470 

HD 57682 

ESPADONS 

May 4'" 2009 

06:05:55 

4x600 

500-1130 


ESPADONS 

May 5 tb 2009 

06:22:17 

4x540 

450-1000 


ESPADONS 

May T* 2009 

06:46:41 

4 x 540 

350-700 


ESPADONS 

May 8"* 2009 

06:21:09 

4x540 

350-800 


ESPADONS 

May 9 th 2009 

06:20:36 

4x540 

300-600 

rSco 

ESPADONS 

May 23 rd 2005 

09:14:25 

4x60 

700-1550 


ESPADONS 

Sep 20 Ul 2005 

05:01:56 

4x30 

600-1200 


ESPADONS 

Aug 5* 2006 

05:25:22 

4x30 

500-1150 


ESPADONS 

Mar 4 ,h 2007 

16:23:15 

4x30 

700-1350 


ESPADONS 

Jun 29 th 2008 

08:03:07 

4x45 

400-1050 



